A Bayesian method for evaluating and discovering disease loci associations

Background: A genome-wide association study (GWAS) typically involves examining representative SNPs in individuals from some population. A GWAS data set can concern a million SNPs and may soon concern billions. Researchers investigate the association of each SNP individually with a disease, and it is becoming increasingly commonplace to also analyze multi-SNP associations. Techniques for handling so many hypotheses include the Bonferroni correction and recently developed Bayesian methods. These methods can encounter problems. Most importantly, they are not applicable to a complex multi-locus hypothesis which has several competing hypotheses rather than only a null hypothesis. A method that computes the posterior probability of complex hypotheses is a pressing need. Methodology/Findings: We introduce the Bayesian network posterior probability (BNPP) method which addresses the difficulties. The method represents the relationship between a disease and SNPs using a directed acyclic graph (DAG) model, and computes the likelihood of such models using a Bayesian network scoring criterion. The posterior probability of a hypothesis is computed based on the likelihoods of all competing hypotheses. The BNPP can not only be used to evaluate a hypothesis that has previously been discovered or suspected, but also to discover new disease loci associations. The results of experiments using simulated and real data sets are presented. Our results concerning simulated data sets indicate that the BNPP exhibits both better evaluation and discovery performance than does a p-value based method. For the real data sets, previous findings in the literature are confirmed and additional findings are found. Conclusions/Significance: We conclude that the BNPP resolves a pressing problem by providing a way to compute the posterior probability of complex multi-locus hypotheses. A researcher can use the BNPP to determine the expected utility of investigating a hypothesis further. Furthermore, we conclude that the BNPP is a promising method for discovering disease loci associations. © 2011 Jiang et al

Expression profiling of familial breast cancers demonstrates higher expression of FGFR2 in BRCA2-associated tumors

BackgroundBRCA1- and BRCA2-associated tumors appear to have distinct molecular signatures. BRCA1-associated tumors are predominantly basal-like cancers, whereas BRCA2-associated tumors have a predominant luminal-like phenotype. These two molecular signatures reflect in part the two cell types found in the terminal duct lobular unit of the breast. To elucidate novel genes involved in these two spectra of breast tumorigenesis we performed global gene expression analysis on breast tumors from germline BRCA1 and BRCA2 mutation carriers. Methodology Breast tumor RNAs from 7 BRCA1 and 6 BRCA2 mutation carriers were profiled using UHN human 19K cDNA microarrays. Supervised univariate analyses were conducted to identify genes differentially expressed between BRCA1 and BRCA2-associated tumors. Selected discriminatory genes were validated using real time reverse transcription polymerase chain reaction in the tumor RNAs, and/or by immunohistochemistry (IHC) or by in situ hybridization (ISH) on tissue microarrays (TMAs) containing an independent set of 58 BRCA1 and 64 BRCA2-associated tumors. Results Genes more highly expressed in BRCA1-associated tumors included stathmin, osteopontin, TGFβ2 and Jagged 1 in addition to genes previously identified as characteristic of basal-like breast cancers. BRCA2-associated cancers were characterized by the higher relative expression of FGF1 and FGFR2. FGFR2 protein was also more highly expressed in BRCA2-associated cancers (P = 0.004). SignificanceBRCA1-associated tumours demonstrated increased expression of component genes of the Notch and TGFβ pathways whereas the higher expression of FGFR2 and FGF1 in BRCA2-associated cancers suggests the existence of an autocrine stimulatory loop

Interplay between estrogen receptor and AKT in estradiol-induced alternative splicing.

BACKGROUND: Alternative splicing is critical for generating complex proteomes in response to extracellular signals. Nuclear receptors including estrogen receptor alpha (ERα) and their ligands promote alternative splicing. The endogenous targets of ERα:estradiol (E2)-mediated alternative splicing and the influence of extracellular kinases that phosphorylate ERα on E2-induced splicing are unknown. METHODS: MCF-7 and its anti-estrogen derivatives were used for the majority of the assays. CD44 mini gene was used to measure the effect of E2 and AKT on alternative splicing. ExonHit array analysis was performed to identify E2 and AKT-regulated endogenous alternatively spliced apoptosis-related genes. Quantitative reverse transcription polymerase chain reaction was performed to verify alternative splicing. ERα binding to alternatively spliced genes was verified by chromatin immunoprecipitation assay. Bromodeoxyuridine incorporation-ELISA and Annexin V labeling assays were done to measure cell proliferation and apoptosis, respectively. RESULTS: We identified the targets of E2-induced alternative splicing and deconstructed some of the mechanisms surrounding E2-induced splicing by combining splice array with ERα cistrome and gene expression array. E2-induced alternatively spliced genes fall into at least two subgroups: coupled to E2-regulated transcription and ERα binding to the gene without an effect on rate of transcription. Further, AKT, which phosphorylates both ERα and splicing factors, influenced ERα:E2 dependent splicing in a gene-specific manner. Genes that are alternatively spliced include FAS/CD95, FGFR2, and AXIN-1. E2 increased the expression of FGFR2 C1 isoform but reduced C3 isoform at mRNA level. E2-induced alternative splicing of FAS and FGFR2 in MCF-7 cells correlated with resistance to FAS activation-induced apoptosis and response to keratinocyte growth factor (KGF), respectively. Resistance of MCF-7 breast cancer cells to the anti-estrogen tamoxifen was associated with ERα-dependent overexpression of FGFR2, whereas resistance to fulvestrant was associated with ERα-dependent isoform switching, which correlated with altered response to KGF. CONCLUSION: E2 may partly alter cellular proteome through alternative splicing uncoupled to its effects on transcription initiation and aberration in E2-induced alternative splicing events may influence response to anti-estrogens.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Low-risk susceptibility alleles in 40 human breast cancer cell lines

Background: Low-risk breast cancer susceptibility alleles or SNPs confer only modest breast cancer risks ranging from just over 1.0 to 1.3 fold. Yet, they are common among most populations and therefore are involved in the development of essentially all breast cancers. The mechanism by which the low-risk SNPs confer breast cancer risks is currently unclear. The breast cancer association consortium BCAC has hypothesized that the low-risk SNPs modulate expression levels of nearby located genes. Methods: Genotypes of five low-risk SNPs were determined for 40 human breast cancer cell lines, by direct sequencing of PCR-amplified genomic templates. We have analyzed expression of the four genes that are located nearby the low-risk SNPs, by using real-time RT-PCR and Human Exon microarrays. Results: The SNP genotypes and additional phenotypic data on the breast cancer cell lines are presented. We did not detect any effect of the SNP genotypes on expression levels of the nearby-located genes MAP3K1, FGFR2, TNRC9 and LSP1. Conclusion: The SNP genotypes provide a base line for functional studies in a well-characterized cohort of 40 human breast cancer cell lines. Our expression analyses suggest that a putative disease mechanism through gene expression modulation is not operative in breast cancer cell lines

Familial hypercholesterolaemia in children and adolescents from 48 countries: a cross-sectional study

Background Approximately 450 000 children are born with familial hypercholesterolaemia worldwide every year, yet only 2·1% of adults with familial hypercholesterolaemia were diagnosed before age 18 years via current diagnostic approaches, which are derived from observations in adults. We aimed to characterise children and adolescents with heterozygous familial hypercholesterolaemia (HeFH) and understand current approaches to the identification and management of familial hypercholesterolaemia to inform future public health strategies. Methods For this cross-sectional study, we assessed children and adolescents younger than 18 years with a clinical or genetic diagnosis of HeFH at the time of entry into the Familial Hypercholesterolaemia Studies Collaboration (FHSC) registry between Oct 1, 2015, and Jan 31, 2021. Data in the registry were collected from 55 regional or national registries in 48 countries. Diagnoses relying on self-reported history of familial hypercholesterolaemia and suspected secondary hypercholesterolaemia were excluded from the registry; people with untreated LDL cholesterol (LDL-C) of at least 13·0 mmol/L were excluded from this study. Data were assessed overall and by WHO region, World Bank country income status, age, diagnostic criteria, and index-case status. The main outcome of this study was to assess current identification and management of children and adolescents with familial hypercholesterolaemia. Findings Of 63 093 individuals in the FHSC registry, 11 848 (18·8%) were children or adolescents younger than 18 years with HeFH and were included in this study; 5756 (50·2%) of 11 476 included individuals were female and 5720 (49·8%) were male. Sex data were missing for 372 (3·1%) of 11 848 individuals. Median age at registry entry was 9·6 years (IQR 5·8–13·2). 10 099 (89·9%) of 11 235 included individuals had a final genetically confirmed diagnosis of familial hypercholesterolaemia and 1136 (10·1%) had a clinical diagnosis. Genetically confirmed diagnosis data or clinical diagnosis data were missing for 613 (5·2%) of 11 848 individuals. Genetic diagnosis was more common in children and adolescents from high-income countries (9427 [92·4%] of 10 202) than in children and adolescents from non-high-income countries (199 [48·0%] of 415). 3414 (31·6%) of 10 804 children or adolescents were index cases. Familial-hypercholesterolaemia-related physical signs, cardiovascular risk factors, and cardiovascular disease were uncommon, but were more common in non-high-income countries. 7557 (72·4%) of 10 428 included children or adolescents were not taking lipid-lowering medication (LLM) and had a median LDL-C of 5·00 mmol/L (IQR 4·05–6·08). Compared with genetic diagnosis, the use of unadapted clinical criteria intended for use in adults and reliant on more extreme phenotypes could result in 50–75% of children and adolescents with familial hypercholesterolaemia not being identified. Interpretation Clinical characteristics observed in adults with familial hypercholesterolaemia are uncommon in children and adolescents with familial hypercholesterolaemia, hence detection in this age group relies on measurement of LDL-C and genetic confirmation. Where genetic testing is unavailable, increased availability and use of LDL-C measurements in the first few years of life could help reduce the current gap between prevalence and detection, enabling increased use of combination LLM to reach recommended LDL-C targets early in life. Funding Pfizer, Amgen, Merck Sharp & Dohme, Sanofi–Aventis, Daiichi Sankyo, and Regeneron

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Gravitational Waves and Gamma-Rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A

On 2017 August 17, the gravitational-wave event GW170817 was observed by the Advanced LIGO and Virgo detectors, and the gamma-ray burst (GRB) GRB 170817A was observed independently by the Fermi Gamma-ray Burst Monitor, and the Anti-Coincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory. The probability of the near-simultaneous temporal and spatial observation of GRB 170817A and GW170817 occurring by chance is $5.0\times {10}^{-8}$. We therefore confirm binary neutron star mergers as a progenitor of short GRBs. The association of GW170817 and GRB 170817A provides new insight into fundamental physics and the origin of short GRBs. We use the observed time delay of $(+1.74\pm 0.05)\,{\rm{s}}$ between GRB 170817A and GW170817 to: (i) constrain the difference between the speed of gravity and the speed of light to be between $-3\times {10}^{-15}$ and $+7\times {10}^{-16}$ times the speed of light, (ii) place new bounds on the violation of Lorentz invariance, (iii) present a new test of the equivalence principle by constraining the Shapiro delay between gravitational and electromagnetic radiation. We also use the time delay to constrain the size and bulk Lorentz factor of the region emitting the gamma-rays. GRB 170817A is the closest short GRB with a known distance, but is between 2 and 6 orders of magnitude less energetic than other bursts with measured redshift. A new generation of gamma-ray detectors, and subthreshold searches in existing detectors, will be essential to detect similar short bursts at greater distances. Finally, we predict a joint detection rate for the Fermi Gamma-ray Burst Monitor and the Advanced LIGO and Virgo detectors of 0.1-1.4 per year during the 2018-2019 observing run and 0.3-1.7 per year at design sensitivity

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO's second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=3.47×10-25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering. © 2019 American Physical Society

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO’s second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h95%0=3.47×10−25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering